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It  was  shown  in  the  preceding  paper ^ that  an  optimum  receiver 
for  a signal  subject  to  temporal  distortions  is  an  adaptive  correlator  which 
maximises  the  quantity 
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where  f . (^wj  = time  t = -yy  of  the  received  waveform  f{t)  s 

where  f(t)  = n(t)  + r s(t)  and 
n(t)  - additive  noise 
r = attenuation  coefficient 
s(t)  = signal  transmitted, 

•ij'j^j  = sample  of  signal  e<t)  at  time  t - 

W « bandwidth  of  receiver  and  s(t) 

To  ~ ti*ne  delay  due  to  round-trip  transit  time  of  signal, 
averaged  over  signal  duration 
*0  = average  time  base  scale  factor  imposed  by  motion 

of  target. 

Ti|zV]  ~ **naple  at  time  t = of  the  time  base  distortion  factor 

T(t)  imposed  by  motion  of  the  medium  and:  target 

accelerations 

c 

Pi  j * cofactor  of  the  (i.j)  component  of  the  normalised 
correlation  matrix  for  v(t)  . 

[p|  = determinant  of  the  normalised  correlation  matrix  for  T(t)  , 
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Obviously  it  is  immaterial  as  to  whether  the  received  waveform  is  adjusted 
to  the  transmitted  signal,  as  stated  here  and  instrumented  in  Fig  3,  or  the 
transmitted  signal  adjusted  to  the  received  waveform,-  as  specified  in  Eq»  (1). 


I Adaptive  Deltic  Correlator 


A.  Conventional  Deltic 

The  usual  Deltic  cross  correlator  ^ uses  two  or  more  recirculating 
delay  lines  for  dynamic  storage  and  multiplication  of  a reference  function 
and  a signal  function  in  a time -compressed  format.  Figure  1 shows  the 
essential  components  of  a Deltic  in  block  diagram  form.  The  upper  delay 
line  of  Fig.  1 receives  hard  clipped  samples  of  tbe  signal  function  at  sample 
intervals  T = (N  + 1)  A , where  N is  the  total  number  of  samples  which 
can  be  stored  on  the  line  at  one  time  and  A is  the  time  interval  per  sample 
on  the  line.  The  output  of  the  line  is  connected  back  to  the  input  so  that  at  a 
time  NA  seconds  later  the  line  will  contain  exactly  the  same  information 
in  the  same  order.  At  a time  (N+  1)  A seconds  later  all  samples  have 
moved  to  the  right  one  sample  position  and  a new  sample  is  received  from 
the  left,  replacing  the  oldest  sample  in  the  line.  Thus  the  signal  function  in 
this  line  processes  by  an  amount  A to  the  right  on  each  complete  circulation 
in  the  sample  interval  T = (N  + 1)  A , and  is  updated  by  receiving  one  new 
sample  on  each  complete  circulation.  The  lower  delay  line  in  Fig.  1 holds 
the  reference  function,  but  this  function  is  neither  processed  nor  updated  as 
it  makes  each  circulation  This  line  is  exactly  (N+  1)  A seconds  long,  so 
that  the  information  merely  circulates  without  precessing  in  the  signal  sample 
interval  T = (N  + 1)  A . The  0 and  1 outputs  o£  the  two  lines  are  connected 
to  two  AND  gates,  which  in  turn  are  connected  through  an  OR  gate  to  a 
counter.  A pulse  input  to  the  counter  signifies  an  identical  polarity  (0  or  1) 

V.  C.  Anderson,  Technical  Memorandum  No.  37  (NR-0 14-903),  Harvard 
University  Acoustics  Research  Laboratory,  January  5,  1956. 
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at  the  instantaneous  output  of  the  two  lines  Tliif  is  a digital  multiplication 


of  the  form  r (nA)  s (nA+  t ) , where  n can  hav©  any  value  from  1 to  (N  -t  1) 
and  represents  the  state  of  relative  precession  between  the  two  delay 
lines.  The  counter  output  provides: 

N>1 

4>(T0)  = /j  r(nA)s(nA  + To)  (2) 

n- 1 

for  each  circulation  of  (N  + 1)A  seconds.  Thus  a correlation  coefficient 
$ (t  ) for  hard-clipped  signals  is  formed  in  each  sample  interval  T = jN*1)A  . 
where  T is  adequately  short  for  proper  sampling  (less  than  one  half  the 
period  of  the  highest  frequency  component  of  the  incoming  signal  function) 

B Adaptive  Deltic 

The  Dsliic  correlator  can  be  modified  to  force  it  to  act  as  an  adaptive 
correlator,  as  shown  in  Fig  2 Here  the  signal  'unction  delay  line  has  » 
number  of  closely  spaced  outputs,  one  for  each  bit  position  near  the  end  of 
the  line.  Adjustment  of  the  signal  function  timebase  is  accomplished  by 
varying  the  switch  position  with  the  proper  rate,  amplitude,  and  polarity 
(with  respect  to  the  center  switch  tap)  to  remove  the  timebase  distortions 
imposed  by  the  medium  as  the  signal  function  iB  reinserted  at  the  beginning 
of  the  delay  line.  Because  the  sampling  interval  T = {N  + 1)A  is  short 
compared  to  the  characteristic  periods  of  the  signal  function  and  the  signal 
samples  are  therefore  redundant,  the  effect  of  the  switch  adding  or  removing 
samples  from  the  signal  will  not  be  detrimental  to  the  correlation  function. 
Improvement  in  the  correlation  output  at  point  A of  Fig.  2 serves  as  a 
monitor  on  the  progress  of  the  correcting  action  of  the  switch.  The  signal 


and  reference  functions  are  recirculated  until  the  output  at  point  A is 
large  enough  to  stop  the  iterative  correction  and  provide  a true  adaptive 
correlator  output. 

The  exact  operation  of  the  Adaptive  Deltic  is  better  described  in 
terms  of  the  more  detailed  diagram  of  Fig.  3.  Each  output  tap  on  the 
signal  function  delay  line  ia  connected  to  a digital  multiplier  of  the  type 
shown  in  Fig.  1,  described  above,  whose  other  input  comes  from  the  ref" 
erence  delay  line  output.  Each  multiplier  feeds  an  integrator  and  envelope 
detector  circuit  of  adjustable  time  constant  T*  . The  outputs  of  all  the 
integrators  and  envelope  detectors  provide  multiple  inputs  to  a "greatest  of" 
selector  of  the  type  shown  in  Fig,  4,  The  most  positive  signal  at  the  control 
grids  of  the  tubes  in  this  circuit  will  asstime  control  by  cutting  off  all  other 
tube sc  and  a coded  output  from  this  circuit  indicates  which  of  the  various 
signal  delay  line  output  taps  provided,  in  conjunction  with  the  reference 
function  output,  the  largest  short  term  correlation  coefficient 

♦ (t  )=  Y/  *(nA)s(nA+T  ) . <3) 

where  N'A  represents  the  adjustable  time  constant  of  the  integrators. 

This  coded  output,  which  measures  the  location  of  the  proper  output  tap  on 
S£  and  therefore  the  timebase  distortion  factor  r(t)  , is  sent  to  a third 
recirculating  delay  line  with  total  delay  equal  to  (N  + 1)A  , This  delay 
line  is  called  the  T(t)  memory,  and  serves  to  program  the  electronic 
switch  S2  on  subsequent  recirculations  so  as  to  control  the  effective  length 
of  the  signal  function  delay  line  in  accordance  with  the  just  measured  values 
of  t (t)  . Since  the  information  in  the  T(t)  memory  circulates  synchronously 


) 
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with  the  signal  function,  its  coded  output  can  operate  switch  S2  in  such  a 
way  that  the  local  timebase  distortions  of  the  signal  function  are  corrected 
instantaneously  (through  the  electronic  switching  operation)  as  they  are  fed 
back  to  the  signal  function  delay  line  input.  Thus  it  is  possible  to  extend  or 
compress  the  effective  length  of  the  signal  function  delay  line  under  control 
of  the  T(t)  memory.  The  correction  is  not  fully  accomplished  in  one  cir- 
culation of  the  lines;  rather,  it  is  an  iterative  operation  in  which  improve- 
ments are  made  on  the  signal  function  on  each  circulation  until  a maximum 
value  is  attained,  and  no  further  improvement  occurs,  in  the  long  term  cor- 
relation coefficient  at  the  output  point  A , Fig  3. 

The  steps  in  the  iterative  process  proceed  as  follows: 

1 Scanning  operation.  The  tap  at  NA  on  the  signal  function  delay 
line  is  connected  to  one  pole  of  switch  for  normal  Deltic  operation  until 
sufficient  precession  of  the  distorted  signal  function  has  occurred  to  bring 
it  into  approximate  alignment  with  the  reference  function. 

2-  Mixed  integration  operation  When  the  sum  #(t  ) of  any  sequential 
combination  of  the  short  term  (T;  = NfA)  integrator  outputs  .),  Eq,  3, 

1 Of  1 

performed  over  the  full  signal  span  by  the  T = {N  + 1)A  integrator, 

<N+1)/NC 

*<To> • <*> 

where  tq  s (tq  = average  over  the  full  signal  time  span  of  the  short 
9 i 

term,  coherent  correlation  time  displacements 


becomes  large  enough  to  exceed  a predetermined  threshold,  even  though 
widely  separated  taps  are  used  for  the  short  term  integrations  «J>.  <t 
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making  up  the  sum,  Sj  is  switched  to  its  other  position  in  which  it  will 
accept  a feedback  from  the  tap  of  S2  that  is  selected  by  the  r(t)  memory. 
This  summation  is  performed  by  the  summing  amplifier  (Fig.  4)  whose  inputs 
are  obtained  one  at  a time  from  the  gates  of  the  "greatest  of"  circuit  Only 
that  gate  connected  to  a plate  circuit  below  ground  potential  will  allow  the 
corresponding  grid  signal  to  pass  to  the  summing  amplifier.  This  incoherent 
summing  of  many  short  integration  time,  coherent  correlations  is  sometimes 
referred  to  as  a mixed  integration.  It  is  used  here  as  a partial  improvement 
on  the  usual  long  term  coherent  correlation  as  a means  of  deciding  that  a 
signal  may  be  present  and  a change  to  the  iterative  procedure  ct  timebase 
adjustment  is  warranted.  The  T integrator.  Fig,  3,  is  followed  by  a dis- 
criminator  which  operates  Sj  when  the  mixed  integration  sum  ®(T<)) 

Eq.  4,  exceeds  the  predetermined  threshold.  By  virtue  of  the  taps  on  S 

2 

this  mixed  integration  has  the  a ^vantage  of  achieving  slightly  different  t 
time  displacements  in  the  sum  of  Eq  4,  and  therefore  one  characteristic 
of  a distorted  time  base.  The  mixed  integration  also  decreases  the  frequency 
resolution  of  the  signal  by  using  much  shorter  integration  times  NA  in  each 
of  the  4>i(Tc  3 ) , Eq  3,  and  therefore  lessens  the  attenuation  of  the  cor- 
relation peak  which  would  otherwise  result  from  small  distortions  «f  the  time 
base.  However,  in  itself  the  mixed  integrator  is  not  an  optimum  receiver,  as 
defined  by  Eq,  1.  because  it  includes  only  displacement  effects  and  not  a 
rubberising,  or  T(t)  , distortion  of  the  timebase.  Furthermore,  it  requires 
no  continuity  of  the  tq  . from  one  short  integration  time  H A to  the  next,  an 
obvious  requirement  for  signal  distortions  imposed  by  the  medium  or  target 
motion.  In  spite  of  these  shortcomings,  it  is  an  excellent  indicator  that  the 
signal  is  so  aligned  with  respect  to  the  reference  that 


a significant  correlation 


coefficient  might  be  achieved  by  subsequent  adjustment  of  the  signal  time- 


3 Timeb&se  adjustment  through  iteration.  With  switch  in  its 
lower  position,  the  timebase  of  the  received  signal  function  is  adjusted 
continuously  by  under  control  of  the  recirculating  T(t)  memory  The 
r{t)  memory  is  constantly  revised  as  the  received  signal  function,  reference 
function,  ?n d r(t)  delay  lines  recirculate  synchronously  and  the  signal  time- 
base  distortions  are  progressively  removed.  On  each  recirculation  only  the 
incremental  change  [ Tk(t)  - Tk  j(t)  ] . in  T(t)  . vdth  respect  to  its  form  on 
the  preceding  (k-1)  circulation,  is  detected  by  the  "greatest  of"  selector, 

Fig.  4.  As  shown  in  Fig,  3,  this  change  is  converted  to  a voltage  level  in  the 
encoder,  to  digital  form  in  the  A/D  converter,  and  then  added  to  ^(t}  , 
the  T(t)  function  on  the  preceding,  (k-1)  , circulation  The  digitised  output 

of  the  adder,  T^t}  , which  is  sent  to  the  T(t)  memory,  is  not  restricted  by  the 
frequency  restraints  of  the  medium  or  the  target  motion  It  is  merely  a 
cumulative  digital  record  of  the  time  varying  outputs  of  the  "greatest  of" 
selector.  However,  the  speed  with  which  electronic  switch  changes 

from  one  tap  to  another  is  set  by  the  r amplifier  (Fig,  3)  frequency  response, 
which  in  turn  is  matched  to  p of  Eq  (1),  The  amplitude  of  the  switching 
motion,  or  the  frequency  with  which  S^,  connects  to  extreme  tap  positions, 

corresponding  to  large  [ Tk(t)  - i^^ft)  1 • is  controlled  by  allowing  the  known 

«r.Z 

— 5-  ratio  to  establish  the  shape  of  the  non-linear,  compressed  (e  g,,  loga- 


rithmic) gain  characteristic  of  the  t amplifier.  The  exact  operating  point 
on  this  gain  characteristic  is  set  by  the  output  from  the  adder  after 

passing  through  a D/A  converter.  The  larger  the  instantaneous  value  of 
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the  smaller  i»  the  differential  gain  of  the  amplifier  The  non-linear,  com- 
pressed shape  of  the  gain  characteristic  and  the  location  of  the  operating 
point  on  this  -haract^ristic,  as  controlled  by  , reflect  the  fact  that 
relatively  large  T(t)  values  are  improbable,  and  are  more  likely  to  be 
present  at  the  “greatest  of1'  selector  output  because  of  adaptation  to  the 
additive  noise  than  to  an  actual  signal  timebase  distortion  In  more 
precise  terms,  making  the  left  side  of  Eq  (1)  a maximum  requires  time  -» 
base  adjustment  of  the  received  signal  function  (to  maximite  the  first  term, 
®*q  [ 1 1 ) tempered  by  rejection  of  relatively  large  and  improbable  T(t)  values 

(to  minimise  the  second  term,  Eq  [l]  ),  The  importance  of  the  second  term 

<r  2 


in  Eq.  (1),  and  therefore  the  degree  of  rejection,  is  set  by  the  ratio 


2 


The  role  of  this  ratio  in  Eq  (1)  is  simulated  in  the  Adaptive  Deltic  by  n 

its  control  of  the  gain  characteristic  of  the  t amplifier  A large  a * . 

characteristic  of  abrupt  target  motions  or  large  boundary  motions  of  the 

medium,  and,  or  small  o ^ , characteristic  of  low  noise  content  m the 

n 

received  waveform,  permit  relatively  frequent  excursion®  of  S,  to  extreme 

2 C 

»T 

tap  positions  On  the  other  hand,  a small  2 ratio  leads  to  a much  higher 

«n 

probability  for  to  contact  a centrally  located  tap  Thus  the  tap  positions 
actually  assumed  by  can  differ  significantly  from  those  called  for  by  the 
T(t)  memory  because  of  the  frequency  and  amplitude  response  of  the  T(t) 
amplifier 

The  actual  correlation  output  is  presented  continuously  at  the  output 
point  A.  Fig,  3.  As  a signal  timebase  adjustment  proceeds  from  one  iteration 
(recirculation  around  the  delay  lines)  to  the  next,  the  output  at  point  A should 
increase  in  magnitude  This  progressive  incre'- se  can  be  used  to  lengthen  the 
time  constant  T"  of  all  the  short  term  integrators  feeding  the  ‘-greatest  of" 


y 


selector  This  is  equivalent  to  lengthening  the  coherent  integration  times  in 


eacl,  of  the  4>.  , Eq.  3,  which  is  justified  as  the  timebase  adjustment  provides 

progressively  longer  coherence  times  of  the  received  waveform  with  respect 
to  the  reference  function.  Eventually  the  short  term  time  constant  T?  should 
become  equal  to  the  long  term  signal  duration  (in  compressed  Deltic  form) 

T , and  switch  should  stay  connected  to  the  central,  (h?  *■  1)A  , feedback 
tap  as  the  signal  timebnse  becomes  fuliy  corrected.  The  average  relative 
velocity  of  the  target,  and  therefore  the  value  of  kQ  in  Eq  (l),  can  be  found 
from  the  average  slope,  over  the  signal  duration  T , of  the  v(t)  function 
stored  in  the  T(t)  memory, 

4 End  of  iteration  The  iteration  process  is  brought  to  a stop  when: 

a)  ceases  to  change  tap  positions,  or 

b)  the  output  at  point  A stops  increasing,  or 

c)  the  output  at  point  A exceeds  a predetermined  threshold. 
All  three  of  these  conditions  satisfied  simultaneously  indicate  a successful 
adaptive  match  of  received  and  reference  functions,  and  therefore  the  presence 
of  a target.  The  following  significant  outputs,  in  addition  to  the  correlation 
coefficient  at  point  A,  are  also  available 

a")  Range  to  target,  as  indicated  by  the  time  at  which 
switched  to  lower  position  to  start  the  iteration  and 
stop  the  regular  precession  of  the  signal  function 

b')  Average  relative  target  velocity,  as  indicated  by  the 
average  slope  of  the  stored  T(t)  function 
However,  it  is  also  possible  that  one  of  the  conditions  (a)  or  (b)  may  occur 
without  a completely  satisfactory  correlation  of  received  and  reference  wave- 
forms This  could  occur  if  distortions  other  than  those  of  the  timebase  such 
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as  signal  fading,  are  present  In  such  an  event,  the  signal  function  may  adjust 
satisfactorily  over  isolated  portions  of  its  timebase  but  these  isolated  portions 
may  not  connect  smoothly  to  provide  a continuous  function  T(t)  In  the  "bad" 


intervals  the  "greatest  of"  selector  either  fails  to  pass  any  signal  exceeding 
the  threshold  of  the  discriminator,  or  those  T(t)  values  passed  are  random 
and  so  widely  separated  that  they  correspond  to  impossible  temporal  changes 


in  the  medium  Then  it  is  desirable  to  construct  a continuous  T{t)  function 
by  interpolating  between  those  portions  of  the  timebase  which  provide  a strong 
and  unambiguous  indication  of  f(t)  in  the  "greatest  of"  selector  To  do  this, 
it  is  necessary  to  know  the  value  of  T(t)  at  the  start  of  the  next  "good"  portion 
of  the  timebase  when  finishing  the  last  ’good"  portion,  the  two  being  separated 
by  a "bad"  interval  This  can  be  done  by  storing  the  T(t)  value  at  the  begin- 
ning of  the  oldest  "good"  portion  of  the  timebase  (that  "good"  portion  which  is 
about  to  appear  next  at  the  output  of  the  T(t)  memory  by  virtue  of  the  recir- 
culation), and  sending  it  to  the  f(t)  amplifier  as  a bias  level  to  be  approached 
during  the  "bad"  interval  The  rate  of  approach,  or  of  switching  S ^ . will  be 
determined  by  the  difference  in  T(t)  at  the  beginning  and  end  of  the  "bad" 
interval  and  by  the  length  of  the  "bad"  interval  This  linear  interpolation  wii* 
allow  to  be  in  a position  to  accept  the  T(t)  values  for  the  next  "good" 
interval,  the  first  part  of  which  would  have  been  lost  to  it  if  it  had  to  approach 
these  well  established  T(t)  values  (S^  tap  positions)  suddenly  through  the  low 
frequency  response  of  the  T(t)  amplifier  Obviously  the  output  at  point  A will 
be  reduced  in  this  case,  and  a decision  must  be  made  with  regard  to  (c)  above, 
the  threshold  level,  as  to  whether  a target  is  present  or  not  However,  the 
adaptive  receiver  is  still  an  optimum  one  for  correcting  timebase  distortions 
whether  or  not  multiplicative  noise  (fading,  etc.  ) is  present  on  the  signal 


Since  the  timebase  adjustment  involves  an  iterative  operation,  in 
general  it  cannot  be  carried  out  in  real  time.  However,  in  many  detection 
situations  sufficient  time  is  available  between  transmissions  to  operate  in 
off  line  fashion  when  examining  probable  target  returns  The  off-line 
operation  is  begun,  and  the  relative  precession  of  the  received  signal  and 
reference  delay  lines  stopped,  only  when  the  mixed  integration  output  of  the 
"greatest  of"  selector  indicates  a target  may  be  present  The  operation  returns 
to  a real  time  search,  and  switch  assumes  its  upward  (original)  position, 
whenever  any  of  the  conditions  (a),  (b)  or  (c)  occur  If  all  three  have  occurred, 
a target  is  declared  present. 

Although  the  T(t)  memory  delay  line  and  the  other  delay  lines  have 
approximately  the  same  length,  the  binary  coding  of  the  T(t)  data  requires 
more  than  one  sample  position  for  each  sample,  whereas  the  hard  clipped 
signal  and  reference  samples  are  fully  described  by  one  bit  each  However, 
this  does  not  require  a higher  packing  density  on  the  r(t)  memory  line 
because  the  highest  frequency  comp  .vue.il  in  T(t)  is  much  lower  than  the 
received  signal  and  reference  upoer  frequencies  Therefore  several  bit 
positions  can  be  assigned  to  a single  T(t)  sample  while  still  maintaining  an 
adequate  sampling  rate. 

It  is  not  intended  that  the  w»y  in  which  o , e * and  o ^ control  the 
1 amplifier.  Fig  3,  should  exactly  simulate  the  role  oi  these  parameters  in 


Eq  (1)  In  general  terms,  p controls  the  time  constant,  or  high  frequency 

2 

, °T 

response,  while  — . controls  the  non  -linear  shape  of  the  gain  characteristic 


of  the  t amplifier,  but  neither  does  so  exactly  in  the  manner  prescribed  by 


Eq  (1)  The  exact  form  of  Eq  (1)  need  not  be  simulated  since  an  iterative 
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procedure  ie  being  followed  to  approach  a maximum  output  level  at  point  A, 

Fig  3,  and  this  maximum  will  occur  for  a particular  final  adjustment  of  the 
signal  timebase  regardless  of  the  detailed  nature  of  the  iterative  steps  which 
approach  it  It  is  only  necessary  that  no  secondary  maxima  occur  at  point  A 
during  the  iterative  approach  to  the  true  maximum,  for  if  such  were  the  case 
the  process  would  stop  once  a secondary  maximum  was  achieved.  However, 
when  a target  is  actually  present  and  the  signal-tc  -noise  ratio  is  adequate  for 
detection  this  is  an  unlikely  event,  for  adjustment  of  the  timebase  of  a target 
return  should  lead  smoothly  to  a true  maximum.  It  is  more  likely  to  occur 
when  the  iterative  process  is  triggered  by  a fortuitous  match  between  receiver 
channel  noise  and  the  reference  function,  in  which  case  the  iteration  will  soon 
stop  by  virtue  of  cause  (a)  or  (b)  of  step  4 above,  and  a target  will  be  declared 
not  present  by  failure  to  exceed  the  threshold,  cause  |c) 
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II  Digital  Computer  Program 


All  of  the  logical  operation*  ascribed  to  the  Adaptive  Deltic  cor- 
relator in  the  previous  section  can  be  programmed  on  a digital  computer 
Moreover,  a general  purpose  computer  allows  a flexibility  of  programming 
which  would  permit  desirable  deviations  from  this  logic  in  specific  cases 
The  computer  has  the  further  advantage  over  the  Adaptive  Deltic  of  treat- 
ing multi  -level  quantising  of  the  received  analog  waveform  fft)  and  the 
reference  function  s(t)  with  greater  ease.  Obviously  the  disadvantage  of 
the  computer  is  its  high  initial  and  operating  costs,  and  for  special  purpose 
operations  such  as  correlation  this  is  seldom  justified  when  the  operations 
are  to  be  repeated  many  times.  Nevertheless,  a digital  computer  with  an 
analog -to -digital  converter  at  the  input  is  a useful  system  for  exploring  the 
capability  of  adaptive  correlation  techniques. 

Operations  which  constitute  major  portions  of  such  a computer 
program  for  adaptive  correlation  are  the  same  as  those  performed  by  the 
Adaptive  Deltic: 

1)  Multiplication  and  integration  of  samples  from  a fixed  and  a pro- 
cessing information  channel,  as  in  a normal  Deltic  program. 

2)  A mixed  integration  consisting  of  the  sum  of  several  outputs  from 
a "greatest  of"  selector  whose  inputs  are  the  envelopes  of  a series 
of  coherent  correlations  involving  integration  times  limited  to  the 
stability  time  of  the  medium  (presumably  shorter  than  the  signal 
duration). 

3)  An  adaptive  change  in  the  timebase  of  f(t)  in  accordance  with  an 
iteration  of  2)  above,  with  the  constraints  of  the  following  quantities 

known  a priori: 
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p , the  normalised  autocorrelation  of  T(t)  , which  determines 
the  permissible  frequency  range  of  t(t) 

b)  <rT2  . the  variance  of  T(t)  . which  determines  the  amplitude 
range  of  T(t)  . 

c)  , the  average  noise  power,  which  determines  the  extent  to 
which  a)  and  b)  above  can  be  allowed  to  constrain  the  timebase 
adjustment  of  f(t)  . 

4)  An  end  of  the  iteration  process  when  either  the  output  correlation 
function  or  the  r(t)  memory  ceases  to  change  significantly. 
Computer  programs  have  already  been  written  for  conventional  Deltic 
correlation,  and  their  extension  to  an  adaptive  operation  is  not  difficult. 


3 C.  N,  Pryor,  "A  High-Speed  Cross  Correlation  and  Fourier  Transform 
Subroutine  Package  for  IBM  7090,"  U S.  Naval  Ordnance  Laboratory 
Memorandum  to -File,  21  June  1962. 
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Ill  Adaptive  Optical  Correlator 


The  same  general  philosophy  but  quite  different  physical  principles  can 
be  used  to  arrive  at  the  conceptual  design  of  the  correlator  which  performs  an 
adaptive  timebase  correlation  by  optical  methods.  Such  an  approach  retains  the 
advantages  of  multi-channel  capability  and  large  time -bandwidth  capacity  char- 
acteristically associated  with  optical  correlators  while  adding  the  adaptive  time- 
base  feature  The  problem  is  somewhat  more  complicated  with  optical  correla- 
tors because  no  two  channels  will  in  general  have  exactly  the  same  T(t)  time- 
base  distortion  factor.  In  fact,  if  the  information  on  the  various  channels  arrives 
at  the  receiver  by  grossly  different  propagation  paths,  or  if  differently  acceler- 
ating targets  appear  on  different  channels,  the  various  T\(t)  , where  the  sub- 
script j refers  to  channel  number,  will  be  completely  independent.  Inauchacase 
it  is  necessary  to  adjust  the  timebases  of  all  channels  independently  and  simul- 
taneously Normally  this  would  seem  to  be  a far  more  difficult  task  than  adjust  - 
ment of  a single  timebase  , as  in  the  Adaptive  Deltic  However,  the  partially 
coherent  panchromatic  optical  correlator  is  particularly  suited  to  this  task  be- 
cause the  time  derivatives  of  the  T.(t)  appear  naturally  at  its  output,  and  these 
can  be  sensed  and  fed  back  to  control  the  writing  of  the  various  signal  functions 
•j\t)  The  partially  coherent  panchromatic  corre  lator  ir  essentially  a combina  - 
tion of  two  optical  correlator  types  described  previously.  1)  the  panchromatic 
optical  correlator*  in  which  time  distortions  in  the  received  signal  are  compen- 
sated by  the  light  source  spectrum,  and  2)  the  partially  coherent  optical  correla- 
tor* in  which  only  partial  light  c oherence  is  maintained  along  the  time  axis  of  the 

4 Ross  E Williams,  "The  Panchromatic  Principle  in  Optical  Filtering" ' Hudson 
Laboratories, Columbia  University  Technical  Memorandum  No  68,  January  30, 
1963) 

4 Ross  E Williams,  ‘’Partially  Coherent  Data  Processing  in  Optical  Systems* 
(Hudson  Laboratories,  Columbia  University  Technical  Memorandum  No  69, 
February  28,  1963) 


receiver  signal  function  so  that  * mixed1  integration  analogous  to  Eq  3 it 
performed  in  the  output  plane. 

Figure  5 shows  a schematic  diagram  of  the  partially  coherent  pan- 
chromatic optical  correlator  Each  point  of  the  area  light  source  illuminates 
only  a restricted  region  of  the  signal  film,  and  the  coherence  of  this  illumina- 
tion falls  off  radially  from  the  center  of  the  region  according  to  a coherence 
factor  y(t-t'1)  which  has  a roughly  Gaussian  shape  (actually  a first  order 
Bessel  function  divided  by  its  argument)  The  totality  of  all  such  overlapping, 
coherently  illuminated  regions  produces  a mixed  integration  at  the  output  of 
the  form 


1 (T0)  = //*  f(t~To)s(t}f(t'-To)s(t?)dtdt‘'  . (5) 

This  result  is  equivalent  to  a superposition  of  many  coherent  correlators, 
each  with  a restricted  aperture  of  tapered  transmitivity  y(t*t  ) overlapping 
one  another,  whose  outputs  are  summed  incoherently.  That  part  of  the  optical 
correlator  to  the  left  of  PI  in  Fig  5 shows  the  partially  coherent  illumina* 
tion,  the  extent  of  whose  coherence  (the  width  of  the  \ (t  t')  factor)  can  be 
varied  by  axial  motion  of  the  lenses  LI  which  causes  each  light  source  point 
in  plane  PO  to  be  imaged  closer  to  or  further  from  plane  PI  , thereby 
changing  the  coherently  illuminated  area  formed  by  the  intersection  of  plane 
PI  with  the  cone  of  light  rays  from  the  source  point 

Thus  the  partial  coherence  principle  produces  a true  mixed  integration 
of  adjustable  coherence  time  which  may  serve,  as  in  the  Adaptive  Deltic  cor- 
relator, as  the  starting  point  for  an  adaptive  technique.  That  part  of  the  optical 

correlator  to  the  right  of  plane  PI  in  Fig  5 is  a conventional  panchromatic 

4 

correlator. 
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A brief  statement  of  some  of  the  principles  of  the  panchromatic 

optical  correlator  is  in  order  to  facilitate  an  understanding  of  its  use  as  an 

adaptive  device.  As  its  name  suggests,  the  panchromatic  correlator  uses  a 

light  source  with  a finite  spectral  width  rather  than  a monochromatic  source 

With  the  received  signal  function  f.(t)  inserted  in  a spatial  plane  and  the 

reference  function  in  a frequency  plane,  the  diffraction  grating  Gj  benes  the 

optic  axis  to  simulate  a frequency  shift  when  the  time  distorted  received 

signals  are  heterodyned  to  a lower  frequency  band  before  recording  on  film 

in  order  to  use  the  low  pass  spectrum  of  the  correlator  more  efficiently. 

The  panchromatic  correlator  produces  a correlation  function  at  its  output 

d-r  (t) 

for  any  T.(t)  distortion  of  the  signals  t.[t~T  (t)]  provided  each  — 

J J J 

is  constant  over  the  optical  aperture  represented  by  y(t-t')  and  regardless 
of  whether  the  T.(t)  are  identical  from  channel  to  channel  In  fact,  each 
frequency  in  the  light  source  correspond*  to  a particular  - . and  thus 
the  various  i\(t)  distortion  factors  will  be  compensated  by  different  pre- 
dominant light  frequencies  in  the  various  channels  at  the  output  These  are 
separated  by  the  color  separation  prism  to  the  left  of  plane  P4  in  Fig  5. 

The  output  format  in  plane  P4  is  an  area  display  as  shown  in  Fig  6,  with 
dr  <t) 

the  distortion  — g* plotted  vertically  for  each  input  channel,  and  range, 

or  real  time,  horizontally  As  the  signal  film  moves,  the  display  changes 
in  real  time,  thereby  providing  the  dependence  (see  Eq  1)  in  the  cor- 
relation function  $(T0)  on  P4  The  record  in  plane  P4  of  a received 
signal  with  variable  time  base  distortion  would  appear  as  a wiggly  trace  of 

greater  intensity  than  the  background,  wandering  among  various  neighboring 
dT.(t) 

— jJ — channels  as  shown  in  Fig.  6,  provided  the  coherence  interval  rep- 


resented by  y(t-t")  , Eq.  5.  is  short  enough  that  the  received  signal  esaen- 


’ 


1 


tially  ha#  a single  timebase  distortion  m each  time  span  Thus,  with  the 

coherence  time  suitably  small,  the  wiggly  trace  conetitutes  a record  of 

instantaneous  timebase  distortions  of  the  received  signal  over  its  time 

duration  A scanning  device  (e  g , a flying  spot  scanner  or  vidicon  camera 

tube)  which  senses  the  shape  of  this  trace  (see  Fig,  7)  converts  the  vertical 
dr.(t) 

or  — deflection  of  the  trace  in  Fig  6 to  a timebase  correction  (speed- 
ing  up  or  slowing  down)  of  the  original  received  signal  function  as  it  is  recorded 
on  film.  A feedback  correction  for  this  purpose  to  a flying  spot  scanner  which 
records  the  received  signal  function  on  film  is  shown  in  Fig  7.  Without  a 
feedback  signal,  the  flying  spot  scanner  provides  an  intensity  modulated  line 
scan  which  is  imaged  across  the  width  of  the  film,  i.e.  , across  the  many 
independent  channels  which  are  to  be  processed  simultaneously.  By  multiplex- 
ing samples  from  each  of  the  independent  receiver  channels  onto  the  control 
grid  of  the  flying  spot  scanner,  a single  scanner  can  expose  all  the  channels 
on  the  film  and  synchronously  demultiplex  them  with  the  cross  channel  line 
scan  shown  in  Fig  7,  The  film  is  advanced  at  a rate  which  causes  each 
successive  line  scan  to  be  imaged  on  the  film  parallel  to  the  last  scan  image 
and  displaced  one  resolvable  distance  along  the  time  axis  from  it 

The  camera  tube  (e  g , vidicon)  scans  the  correlator  output  plane  in 
dT.(t) 

the  vertical  or  — — direction  synchronously  with  the  vertical  line  scan  of 
the  flying  spot  scanner  exposing  the  Rignal  film.  (Alternatively,  a flying  spot 
scanner -photomultiple r combination  can  be  used  to  "read''  the  correlator  out- 
put plane  information  if  it  is  already  recorded  on  film. ) When  the  output  plane 


flying  spot  scanner  exposing  the  signal  film  However,  when  a correlation 

dr.(t) 

peak  occurs  for  a non  rero  — — . the  feedback  signal  to  the  deflection 

plates  of  the  flying  spot  scanner  causes  the  spot  image  to  deflect  along  the 

time  dimension  ilengthwise  dimension;  of  the  Film,  thereby  adjusting  the 

timebase  of  that  particular  channel  to  compensate  ior  the  local  timebase 

distortion  Thus  the  line  scan  is  no  longer  » straight  line  scan  across  the 

width  of  the  film,  but  departs  in  the  time  direction  whenever  necessary  to 

correct  for  recognised  time  distortions.  The  direction  and  magnitude  of  the 

timewise  deflection  of  the  spot  are,  of  course,  set  by  the  sign  and  magnitude 
dT  (t) 

of  the  local  — -g — encountered  in  the  correlator  output  plane,  and  because 
of  the  low  inertia  of  the  CRT  electron  beam  and  high  frequency  response  of 
the  deflection  amplifier  all  signal  channels  can  be  adjusted  simultaneously 
and  independently  As  in  the  Adaptive  Deltic  correlator,  care  must  be  taken 
to  avoid  adapting  the  received  signal  channel  timebases  to  noise  in  the  channels, 
and  the  proper  restraints  incorporated  in  the  low  pass  filter  and  non  linear 
amplifier  in  Fig.  7 serve  the  same  purpose  as  the  equivalent  ones  in  the 
Adaptive  Deitic  schematic,  Fig,  3 

When  the  feedback  technique  is  used  it  is  obvious  that  real-time 
operation  is  no  longer  possible,  for  the  feedback  operates  upon  a pre-existing 
signal  to  correct  it  for  a more  coherent  correlation:  As  the  signal  timebaee 
becomes  better  adapted  to  remove  unwanted  distortions  by  successive  iterations 
through  the  feedback  path,  as  described  above,  the  coherence  interval  |t  t j 
can  be  made  longer  and  longer  until  eventually,  when  all  distortions  have  been 
removed,  it  occupies  the  full  aperture  of  the  correlator-  This  variation  of 
the  coherence  factor  y(t-tp)  is  produced  by  axial  adjustment  of  lenses  LI 
in  Fig-  7 Starting  with  LI  so  positioned  that  the  coherence  area  on  PI  is 
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small  (plane  PO  nearly  imaged  upon  Pi)  the  lenses  LI  are  gradually 
moved  toward  PO  until  they  are  a focal  length  f from  it,  when  parallel 
light  rays  from  any  given  point  in  PO  illuminate  the  full  aperture  of  PI 
coherently  and  a completely  coherent  correlation  is  performed  At  the 
same  time,  the  wiggly  trace  in  the  output  plane  P4  , Fig  6,  should 
straighten,  becoming  completely  straight  when  completely  coherent  cor 
relation  is  achieved. 

The  process  actually  does  proceed  in  real  time  until  a partially 
coherent  correlation  exceeds  the  threshold  level  for  signals  in  plane  P 4 , 
since  no  feedback  ia  present  until  that  time.  When  the  threshold  level  is 
exceeded  and  feedback  occurs,  the  signal  inputs  to  the  flyjng  spot  scanner 
control  grid  must  be  switched  from  the  signal  channels,  where  new  signals 
continue  to  arrive,  to  a recorded  version  of  the  previous  signals  which  first 
gave  a partially  coherent  correlation  indication.  This  is  done  most  readily 
by  recording  all  received  signals  on  a loop  of  magnetic  tape  whose  length 
is  matched  to  the  maximum  time  -aperture  of  the  correlator,  and,  when  the 
threshold  is  exceeded,  replaying  the  loop  into  the  flying  spot  scanner  rather 
than  continuously  erasing  and  recording  new  cignals,  as  in  the  case  with  no 
correlation  present  above  the  threshold 


Tfee  purpose  of  each  of  the  adaptive  correlators  described  above  is 
to  determine  whether  the  relatively  low  level  correlation  resulting  from  a 
partially  coherent  integration  of  the  signal  timebase  actually  is  associated 
with  a target  or  not  If  it  is,  the  correlation  peak  should  ri«e  as  the  time- 
base  i*  adjusted  tor  local  distortions  reaching  a maximum  when  the  cor 
relation  becomes  completely  coherent  Or  the  other  hand,  when  the  par 
tially  coherent  result  exceeds  vhe  first  threshold  only  because  of  a fortuitous 
matching  of  the  received  signal  channel  noise  with  the  reference  function 
(a  not  unlikely  event)  no  amount  of  timebase  adjustment  through  successive 
iterations  will  enhance  the  correlation  peak  when  the  feedback  is  made  wit  i 
proper  frequency  and.  amplitude  restraints  In  this  case  the  iterations  car 
be  cut  off,  and  normal  real-  time  operations  restored,  when  the  correlation 
peak  fails  to  exceed  a second,  higher  threshold  after  a fixed  number  of 
iterations  (step  4b  above)- 

It  is  important  to  note  that  the  problem  in  adaptive  timebsse  correlation 
is  not  primarily  one  of  enhancing  the  signal -to  noise  ratio,  but  rather  one  of 
increasing  the  time  coherence  oi  the  signal  so  that  an  already  adequate  signal 
to  -noise  ratio  will  provide  a suitably  strong  correlation  peak  after  timebai  e 
compensations  have  been  made  Thus  the  technique  is  not  primarily  designed 
for  weak  signal  detections,  but  for  greater  probability  of  defection  when  th  j 
signal  to  noiee  ratio  is  adequate  but  either  target  or  boundary  motions  linr  it 
the  coherence  of  the  signal 
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